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ABSTRACT 


We have revisited the problem of propagation of toroidal and linear Alfven waves 
formulated by Heinemann and Olbert [1980] to compare WKB and non-WKB waves and 
their effects on the solar wind. They considered two solar wind models and showed that 
reflection is important for Alfven waves with periods of the order of one day and longer, and 
that non-WKB Alfven waves are no more effective in accelerating the solar wind than WKB 
waves. There are several recently published papers which seem to indicate that Alfven 
waves with periods of the order of several minutes should be treated as non-WKB waves 
and that these non-WKB waves exert a stronger acceleration force than WKB waves. The 
purpose of this paper is to study the origin of these discrepancies by performing parametric 
studies of the behavior of the waves under a variety of different conditions. In addition, we 
want to investigate two problems that have not been addressed by Heinemann and Olbert, 
namely, calculate the efficiency of Alfven wave reflection by using the reflection coefficient 
and identify the region of strongest wave reflection in different wind models. To achieve 
these goals, we investigated the influence of temperature, electron density distribution, 
wind velocity and magnetic field strength on the waves. The obtained results clearly 
demonstrate that Alfven wave reflection is strongly model dependent and that the strongest 
reflection can be expected in models with the base temperatures higher than 10* K and 
with the base densities lower than 7 x 10 7 cm" 3 . In these models as well as in the models 
with lower temperatures and higher densities, Alfven waves with periods as short as several 
minutes have negligible reflection so can be treated as WKB waves; however, for Alfven 
waves with periods of the order of one hour or longer reflection is significant, requiring a 
non-WKB treatment. We also show that non-WKB, Unear Alfven waves are always less 
effective in accelerating the plasma than WKB Alfven waves. Finally, it is evident from 
our results that the region of strongest wave reflection is usually located at the base of 
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the models, and hence that interpretation of wave reflection based solely on the reflection 
coefficient can be misleading. 


3 



INTRODUCTION 


The existence and importance of Alfven waves in the solar wind has been known and 
extensively studied since the first work on this subject by Belcher and Davis [1971] and 
Belcher [1971]. The main purpose of those studies was to understand the role played by 
Alfven waves in the wind acceleration and their possible association with plasma fluc- 
tuations observed in the solar wind. To gain physical insight into these problems, one 
must consider the propagation of Alfven waves in a highly inhomogeneous medium with a 
background flow. In general, the problem is difficult because Alfven waves propagating in 
inhomogeneous media may be coupled to other magnetohydrodynamic (MHD) waves, and 
then rigorous treatment would require finding solutions to the full set of MHD equations. 
However, early papers on this subject by Ferraro [1954] and Ferraro and Plumpton [1958] 
have shown that under some simplifications, linear Alfven waves can be fully separated 
from other MHD waves. An additional simplification that has been often used to describe 
the propagation of Alfven waves in the solar wind is the WKB approximation, which is an 
appropriate procedure for waves with wavelengths shorter than local characteristic scales 
[e.g., Belcher 1971; Hollweg 1973; Belcher and Olbert 1975; and Jacques 1977, 1978). The 
WKB approximation, however, has long been known to be incorrect for wavelengths of 
typical Alfven waves observed in the solar wind. Therefore, Heinemann and Olbert [1980, 
hereafter H&O] proposed a new mathematical approach which is mdependent of the WKB 
assumption and is based on the full Alfven wave equations derived for the wave velocity 
and magnetic field perturbations [see also Hollweg 1990, and Barnes 1992]. The idea that 
Alfven waves propagating in solar and stellar atmospheres must often be treated as non- 
WKB waves has been recently confirmed by An et al. [1989, 1990], Barkhudarov [1991], 
Velli, Grappin and Mangeney [1991], Velli [1993] and Lou and Rosner [1993]. In addition, 
Moore et al. [1991] have used non-WKB Alfven waves to explain heating in solar coronal 

holes. 
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The fact that Alfven waves propagating in nonuniform atmospheres may become non- 
WKB waves and suffer strong reflection has been known for a number of years [e.g., 
Thomas, 1983; Campos, 1987; and references therein]. Several different methods have been 
developed to address this problem. First, some attempts have been made to obtain full 
analytical solutions to the derived Alfven wave equations and to determine the height in the 
atmosphere where the reflection becomes dominant in the wave behavior. Obviously, the 
closed form solutions can only be found for very few special cases; one celebrated example 
is the propagation of Alfven waves in a plane parallel isothermal atmosphere permeated 
by a uniform and vertical magnetic field [Ferraro and Plumpton, 1958]. In this model, the 
upward propagating Alfven waves of all frequencies are reflected and they always become 
standing waves. Second, it has been suggested that in order to determine the region in the 
atmosphere where wave reflection is the most important, it is sufficient to calculate a local 
cutoff frequency, which is usually defined as a ratio of wave velocity to a characteristic 
scale height. Having obtained the local cutoff frequency, one may then find the height in 
the atmosphere at which the wave frequency is comparable to the cutoff [Rosner, Low and 
Holzer, 1986; Musielak, Fontenla and Moore, 1992] and this height will identify the region 
of strong wave reflection. Finally, numerical solutions to the Alfven wave equation can be 
obtained and non-WKB effects can be studied in detail [see H&O, also Hollweg 1990, and 
references therein]. We shall use this approach because the method presented by H&O is 
both a comprehensive and elegant treatment of Alfven wave propagation in the solar wind. 

There are two main conclusions in H&O’s paper: (1) Alfven waves in the solar wind 
with periods of one day and more must be treated as non-WKB waves, but shorter period 
waves can be described as WKB waves; (2) Non-WKB Alfven waves are no more effective 
in accelerating the solar wind than WKB Alfven waves. The problem is that results 
recently obtained by Barkhudarov [1991], Moore et al. [ 1991 ], Velli et al. [ 1991 ] and 
Velli [1993] seem to be in contradiction with the H&O conclusions. In particular, the 
authors of the three latter papers clearly indicate that Alfven waves with periods as short 
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aa 5 - 20 minutes can be reflected in the solar corona and provide acceleration and heating 
in the solar wind. In addition, Barkhudarov and Velli claim to have shown that some non- 
WKB Alfven waves are much more effective in accelerating the wind than WKB waves. 
It must be noted that most of these recent results have been obtained by using methods 
which differ from the H&O approach; i.e., Barkhudarov and Velli modified the basic set of 
H&O equations by introducing new sets of variables and using different renormalization 
factors. More importantly, slightly different coronal and wind model have been used for 
each calculation. Because H&O performed their calculations for only two solar wind models 
and because Alfven wave reflection depends sensitively on the choice of parameters that 
specify the model, our primary goal in this paper is to extend the H&O results over much 
broader (but observationally plausible) ranges of coronal and solar wind parameters. In 
addition, we want to study the spatial distribution of wave reflection to find the location 
of region in the model where reflection is strongest and to calculate the global efficiency 
of wave reflection by using the reflection coefficient. 

To achieve these goads, we have adopted the H&O approach and numerically solved 
their wave equations. Our first task was to confirm their analysis and results, which we have 
successfully done. Then we investigated the propagation and reflection of Alfven waves 
in a broad spectrum of wind models constructed for different magnetic field strengths, 
temperatures, and densities. The necessity for such investigations becomes apparent from 
a superficial study of H&O. A clear statement of how reflection depends on the complete 
set of coronal parameters used in the H&O model would be an important addition to 
the understanding of Alfven wave reflection in the solar corona. It is also of interest to 
calculate the dependence of the global reflection coefficient on wave frequency in different 
models and compare the results to those previously obtained by other authors. We have 
also found where in the models wave reflection is strongest and compared the acceleration 
force for WKB and non- WKB waves. We believe that our results are helpful in clarifying 
the role played by WKB and non- WKB Alfven waves in the solar wind. 
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BASIC FORMULATION AND PREVIOUS RESULTS 
In general, the problem of propagation of Alfven waves in the solar wind is complicated 
because the waves are usually coupled to other MHD modes. However, as shown by H&O, 
the coupling can be neglected when linear, toroidal Alfven waves are considered. These 
waves have only <f> components in spherical polar coordinates and, according to H&O, they 
can be described by two new wave variables: u = V* + V Ai p and v - V* - V A<t> , where 
and V A <t, represent the ^-component of the wave velocity and the Alfven wave velocity, 
respectively. A model of steady state and axisymmetric solar wind flow is assumed for the 
background medium, with the poloidal Alfven velocity being related to the poloidal flow 

velocity by 

V A = +r? 1/2 V , W 


where rj is the ratio of local gas density p to the gas density />« at the Alfvenic critical point, 
which is defined as a point in the solar wind at which V = V A . Note that the Alfvenic 


critical density p a is a given field line constant. Because the flow is aligned with the local 
magnetic field, H&O introduced the arc length, s, along a given field line and defined two 

new wave variables 


/ = 


1 - 1 


1/2 


ll/4 


U 


9 = 




1/2 


7 1 / 4 


■v , 


(2a) 

( 26 ) 


which represent the inward (/) and outward ($) propagating waves. By using the above 
definitions and these approximations to the set of ideal MHD equations, H&O derived the 
following equations: 




[£ + (K + v,,A].-/(v + i««. 


(за) 

( зб ) 


where * = - ln( ify 1/4 ) with R = r sin 9 being a distance to the symmetry axis. Equations 
(3a, b) showed that the wave amplitudes / and g are coupled together through the density 
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gradient difr/ds and that the propagation of outward (inward) waves is modified by the 
existence of inward (outward) waves; the latter results from wave reflection due to the 
local density gradient. Note that in the WKB approximation, gradients in the background 
medium do not affect the wave propagation and, therefore, the inward propagating waves 
are not present. 

To reduce the set of derived partial differential Equations (3a, b) to ordinary differ- 
ential equations, H&O assumed that both / and g are periodic functions of time and used 
Fourier transforms 


/(M) = f{*)e ,u,t i <?(M) =9{*)e > 


(4) 


where f(s) and <)(.») are the Fourier amplitudes. Using these transforms and introducing 
some new quantities defined below, they showed that Equations (3a, b) can be written in 
the following dimensionless form: 


/ \ (ix (i 


/ \ dx d 


j ,dzd1> 

f = 9(v - u A )-- , 

- ft j_ \ dx dt(> 


where x = r/r a , and 


ds 

dl= — = 

r a 


1 + 




(5a) 


(56) 


( 6 ) 


is the arc length, and 0 = 0(x) represents the equation of the field line. In addition, the 
parameter a is defined by 


a = 


u >r a 2irr a 


rV a ' 


(7) 


with r„ describing the location of the Alfvenic critical point, u and ua bring the flow and 
Alfven velocities normalized by the Alfvenic critical velocity V a , and r being the period of 
Alfven waves, r = 2tt/u>. Because the solutions to Eqs. (5a and b) are complex, we must 
consider |/| and | cj| to get physically meaningful quantities. It is also important to realize 




OH POiftr^UAUSY 


8 



that both |/| and |$| are global quantities such that they are determined at a given height 
by the reflection process taking place in other parts of the flow. We shall discuss this issue 
in more detail in the next section of this paper. 

To solve Equations (5a, b), the boundary conditions have to be specified. One bound- 
ary condition is determined by a singularity in Equation (5a) which occurs at the Alfvemc 
critical point. To remove this singularity, H&O proposed to take / = 0 at the Alfvenic 
critical point as the first boundary condition; this means that there is no wave reflection at 
the critical point and that Alfven waves of all frequencies are moving outward in the Sun s 
frame of reference. Above the critical point, Alfven waves of some frequencies may, however 
still be non-WKB waves. H&O also analyzed in detail the wave action for Alfven waves 
and found that the Fourier amplitudes f(s) and g(s) must obey the following conservation 

relation 

|sf| 2 ~ |/| 2 = con st . (®) 

Combining the first boundary condition with this relation, one obtains $(s) = const at the 
Alfvenic critical point. As a matter of convenience, one may take <j(s) = 1 at the Alfvenic 
critical point as the second boundary condition. It is interesting to note that the boundary 
conditions f{s) = 0 and |(s) = 1 are also WKB solutions to Equations (5a, b). Because 
the differential Equations (5a, b) and the chosen boundary conditions contain only one 
parameter a, the character of the solutions will be determined by this single parameter. 
However, it is important to realize that a depends on the wave frequency, the location of 
the Alfvenic critical point, and the Alfven velocity there, which in turn are determined 
by the magnetic field strength and the electron density distribution. By changing these 
parameters we can investigate the behavior of Alfven waves of different periods propagating 
in different wind models. It should be mentioned that the set of Equations (5a, b) is very 
stiff for large values of a. In some cases it was impossible to obtain results for periods 
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shorter than a few minutes, even with the use of a subroutine especially developed to treat 
stiff equations. 

Two models of the solar wind were used by H&O, namely, the so-called Munro- Jackson 
model [Munro and Jackson 1977], for a polar and nonspherical coronal hole and the “spher- 
ical” model in which the parameters describing the background flow are based on an empir- 
ical density profile. To reproduce the results of H&O, we numerically integrated Equations 
(5a, b) subject to the boundary conditions / = 0 and g = 1 at the Alfvenic critical point. 
We have reconstructed their spherical solar wind model by using the same data adopted 
from Sittler [1978]. The solar wind base was set at r 0 = 1.03.R© and physical parameters at 
r 0 had the following values: magnetic field Bo = 2 G, electron density n eo = 4 x 10 8 cm -3 , 
temperature T 0 = 1.3 x 10 6 K. Then, we calculated the Alfven velocity with distance from 
the Sun and the resulting wind velocity using Eq. (1), with r a = 20 f?©, as suggested 
by H&O. Finally, we have reproduced quantitatively the behavior of Alfven waves in the 
spherical solar wind model by changing the wave period in the range considered by H&O. 
We have found full agreement with their results and were able to draw the following con- 
clusions: (i) all Alfven waves with periods of one day and longer are non-WKB waves 
in the solar wind; (ii) there is an enhancement of the wave energy density for non-WKB 
waves when compared to that found for WKB waves; (iii) there is no corresponding en- 
hancement of the wave energy flux density for non-WKB Alfven waves; and (iv) non-WKB 
Alfven waves are indeed no more effective in accelerating the solar wind than W^KB Alfven 
waves. To test the validity of these conclusions for different solar wind models, we have 
performed a number of parametric studies described in the following two sections. 

WAVE REFLECTION 

To present the results of our parametric studies, we must first describe the wind models 
used in our calculations. We have adapted the method developed by Yeh [1970] to construct 
a series of different wind models. In this method the models are calculated by specifying the 
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physical parameters , i.e. temperature and electron density at the sonic critical point. The 
mass loss rate, M, and the polytropic index, a„, must also be given. All models considered 
in this paper are calculated by taking M = 1.44 x 10 12 <js 1 = 2.3 x 10 M^yr and 
ctp = 1.05; these two parameters remain fixed at these values in our parametric studies. 

For different physical parameters at the sonic critical point, the wind models have different 
distributions of the wind velocity, temperature and density, and also have different values 
of these parameters at the base of the wind which is fixed at r = 1.037*0. This allows 
us to investigate Alfven wave reflection in models that have both the distribution of basic 
physical parameters and their values at the base of the wind significantly different than used 
by H&O. The decision to perform parametric studies in different wind models is justified 
because there are no commonly accepted coronal and wind models based on observational 

(l&t 8»a 

The method applied in this paper to construct the coronal and wind models allows 
us to perform three types of parametric studies. First, we may investigate Alfven wave 
reflection in models characterized by the same density and temperature at the base of the 
wind but having a different value of the magnetic field strength there. Second, we may 
study the dependence of wave reflection on changes in the density at the base by assuming 
that both the base magnetic field and the base temperature are fixed; it must be noted 
that the method chosen here to construct the wind models does not formally allow us to 
fix the temperature at the base and at the same time to vary the base density [Yeh 1970]. 
Therefore our study is restricted to a narrow range of temperatures that changes the base 
density by approximately a factor of 2. Finally, we fixed the magnetic field strength at 
the base and calculated wind models for four different values of the base temperature; as 
mentioned above, this automaticaUy leads to four different values of the density at the 
base. For each of these cases, we solve the wave equations (5a, b) to find the functions 
|/| 2 and \g\ 2 . In addition, to compare our results with those previously obtained, we 
calculate the reflection coefficient. As discussed by Leroy [1980], some caution is required 
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in properly defining the reflection coefficient. Namely, one must clearly identify a region in 
the wind model where wave reflection does not occur, so that Alfven waves existing in this 
region are outward propagating (WKB) waves in the Sun’s frame of reference. Assuming 
that the region exists in the wind model, then the transmission coefficient is defined as 
the ratio of the wave energy carried by the outward propagating waves at the base to that 
carried by these waves in the region without reflection. In the approach considered m this 
paper, Alfven waves of all frequencies have / = 0 at the Alfvenic critical point; beyond 
that point the solar wind convects inward waves outward in the Sun’s frame of reference. 
There may still be reflection beyond the critical point, so that some waves move inward in 
the solar wind frame, but this does not affect the reflection coefficient. Our calculations 
demonstrate that in most of our models waves with periods of the order of one day and 
shorter are also WKB waves beyond the critical point. Therefore, we use the location of 
the Alfvenic critical point to identify the region beyond which all waves propagate outward 
for our definition of the reflection coefficient. To derive the expression for the reflection 
coefficient, we begin with the conservation of wave action, N, written m the following form: 


dN 

+ V • (V g iV) = 0 , 


( 9 ) 


where 

n - pV (-£- ) , (io) 

N - pV \v + v A V- V A ) 

and V g is the group velocity of Alfven waves. After taking the time average of Eq. (9), 
we find that V • (V g JV) = 0 because the time average of the first term vanishes; this can 
be also written as r 2 V g N = const , where N is the averaged wave action. Now, we may 
define the reflection coefficient, 7 Z, for the outward propagating Alfven waves as H 1 T, 
where the transmission coefficient, T , is given by 


T = 


(r 2 V g iVou t )c 
(r 2 V g N 0 Ut)o ' 


(ID 
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with the subscript “c” and “o” denoting the values of the considered quantities at the 
Alfvenic critical point and at the wind base, respectively. Substituting Eq. (10) for the 
outward propagating waves only into Eq. (11) and using V g = V + V A , we find that 
j- _ \/\g 0 1 2 and that 1Z = |/<?| 2 /|y<>| 2 - The reflection coefficient allows estimating the total 
amount of wave energy reflected between the base of the wind model and the location of 
the Alfvenic critical point. 

Before we present the results of our studies, we first want to briefly discuss a simple 
physical picture of wave reflection which emerges from geometrical optics and can be 
helpful in understanding the wave behavior in different wind models. In this picture, a 
nominal wavelength A for an Alfven wave of a given period r is defined at each height in 
the wind model as A = V a t. Because the background medium is inhomogeneous, both 
V A and A must be functions of height. Now, as known from geometrical optics the wave 
propagation is affected by the inhomogeneity of the medium only when the wavelength 
becomes comparable or longer than the characteristic wave velocity scale height. Applying 
this to our case, we have A = H A = V A l{\dV A j dx\), where H A is the Alfven velocity scale 
height; note that H A accounts directly for the density and magnetic field gradients and 
that both the temperature and wind velocity gradients have only indirect effects. It can 
be shown that for any wave-bearing medium, reflection is weak when A << H A and ^at 
it becomes strong when A H A . Obviously, the above criteria are too simple to quantify 
the process of Alfven wave reflection in a wind model in any detail but they are helpful in 
understanding the results based on the solutions of the Alfven wave equations (5a and b) 
and interpreting the behavior of the wave amplitudes |/| and |<j| . 

The results of our parametric studies are presented in six-panel figures (see Figures 1, 
2 and 3). The first three panels in each figure show distributions of density, temperature, 
and Alfven and wind velocity with distance from the Sun. In the next two panels, we 
present changes in the wave (Fourier) amplitude \f\ 2 caused by different wind models and 
wave periods; note that once |/| 2 is known, |j| 2 can always be calculated by using Eq. 
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(8) with const = 1.0. Finally, the last panel shows the reflection coefficient as a functwn 
of wave period. In the following, we discuss in detail the results presented in Figures 1, 2 

and 3. 

In our first parametric study we calculated the temperature and density distributions 
according to the method described above, getting the temperature and the electron density 
at the base to be To = 1.6 x 10* K and «.. = 7.0 x 10 7 cm" 3 , respectively. This we call 
Model 1 (see Figs, la and b). Because the strength of magnetic fields in solar coronal 
holes is presently not well known [e.g., Parker, 1991], we have assumed that the field can 
range from 2 G to 15 G, and preferentially is of the order of 10 G (e.g., Moore et al. 1991]. 
Any change, in the strength of the magnetic field at the solar wind base give different 
Alfven velocities which lead to different locations of the Alfvenic critical point (see Fig. 
lc) and give different models. If the wave frequency u is held constant for all models 
with different magnetic field strengths, then the parameter a is different for each location 
of the Alfvenic critical point and we expect different wave behavior in each model. The 
waves propagating in the model with the strongest magnetic field will have the longest 
wavelengths and, according to the above criteria, will be reflected more than waves of the 
same frequency propagating in other models with weaker magnetic fields. In addition, it 
must be noted that reflection decreases with distance from the model base and becomes 
zero at the point where the Alfven velocity has a maximum; as this point is approached, 
H a - oo, and the medium can be treated locally as homogeneous. Bearing this result in 
mind, we now show the dependence of the wave amplitude \f\ 2 solely on the magnetic field 
strength. We fixed the wave period to be 0.04 d * l k and plotted |/|* versus the distance 
from the base (see Fig. Id). The figure clearly shows that the value of |/| increases with 
increasing magnetic field strength and that it is greatest at the base of each wind model. 
As seen in Figure le, the wave amplitude |/f goes through a local maximum where the 
Alfven velocity peaks (this occurs at the point where d^/dx = 0, see also Eqs. 5a and 
b) and that it is practically zero above 1O* 0 . To understand the result, presented in 
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Figure Id, we must now recall our interpretation of the function |/[ given in the previous 
section. As discussed before, the wave amplitude |/| is not a local quantity but instead 
it is determined by the reflection process taking place everywhere in the model; more 
specifically, because of the boundary condition imposed on / at the Alfvenic critical point, 
each model can be formally divided into the lower (from to the base to the critical point) 
and upper (from the critical point to infinity) part, and the local value of |/| is determined 
by the global wave behavior in each part separately. In other words, the wave amplitude 
|/| 2 reveals how important non-WKB effects are globally in the model (1/1=0 for WKB 
waves) and its value at any height provides information about the global reflection process 
occurring above this height; the local value of |/| 2 does not indicate, however, whether 
reflection is locally strong or weak (see below for more details). Therefore, the fact that 
|/| 2 has a maximum at the same location as the maximum of the Alfven velocity does not 
mean that the maximum of reflection occurs there. It rather indicates that there is a region 
of strong reflection below and above the location of the maximum of the Alfven velocity, 
and that this leads to the formation of the region where waves can be trapped thereby 
causing the calculated enhancement of |/| 2 . To identify the height in the model where 
the strongest wave reflection occurs, one must calculate the derivative of |/| with respect 
to x, take its absolute value, and then plot it versus distance. The maxima and zeros of 
the absolute value of d|/| 2 /dx are located at these heights where wave reflection is locally 
the strongest and weakest, respectively. For example, by taking the real and imaginary 
part of the function |/| 2 and using Eqs. (5a and b), one may find that d\f\ / dz = 0 
at the point where dxfx/dx = 0 (the maximum of the Alfven velocity). This is consistent 
with the physical picture discussed above which requires no reflection at the point where 
dV A /dx = 0. Even though, we do not show the plots of the derivative with distance, the 
results presented in Figure Id clearly indicate that the region of strongest wave reflection 
is located at the base in all considered models. As expected, the weakest wave reflection 
(d|/|/dx = 0) always occurs at the maximum of the Alfven velocity. 


15 



We now investigate the dependence of the wave amplitude |/| J on the wave period 
(see Fig. le). The presented results demonstrate that |/| 2 increases with increasing wave 
periods and for sufficiently long periods reaches a maximum at the peak of the Alfven 
velocity. In the light of our discussion of |/| 2 , this coincidence simply indicates that locally 
wave reflection decreases toward the peak and that it reaches zero at the location of the 
maximum of the Alfven velocity. Recalling now our interpretation of the wave amplitude 
|/| 2 , it is seen from Figure le that long period Alfven (say, 5 days) waves are more strongly 
reflected in the model than those with shorter periods. Note that waves with periods of 5 
days and longer are non-WKB waves even beyond the critical point and that this agrees 
with the discussion in H&O of the behavior of these waves when the parameter a reaches 
a certain value, approximately 0.5 and less. Finally, results presented in Figure If show 
that the reflection coefficient 1Z, defined according to our Eq. (11), also increases with 
increasing magnetic field strengths and wave periods. Even Alfven waves with periods of 
the order of several days are only marginally reflected in our wind model providing that 
the magnetic field is not very strong, ( B a = 2 G). For stronger fields there is always a range 
of wave periods for which the reflection coefficient remains constant. This seems to be in 
contradiction with the results presented in Figures Id and le, which instead demonstrate 
that non-WKB effects should be more prominent when the magnetic field strength and/or 
the wave period are increased. The discrepancy is discussed in more detail in the next 
paragraph of this section. 

The results presented in Figure 2 show the dependence of Alfven wave reflection on 
two wind models (Models 2 and 3, see Figs. 2a and b) which have similar distributions of 
temperature and density but two different values of these parameters at the base. T 0 — 
1.28 x 10 # and 1.35 x 10 6 K\ these two temperatures correspond to the following values 
of the base density: n eo = 7.0 x 10 8 and 3.75 x 10 8 cm' 1 . Note that the distribution 
of temperature and density in both models is very similar and that the main difference 
between the models is the different value of the density at the base; the latter obviously 
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affects the value of the Alfven velocity at the base and leads to slightly larger Alfven 
velocity scale height, H A , in Model 3. This explains the results presented in Figure 2d, 
which show that waves of the same period are more non-WKB waves in Model 2 than in 
Model 3 because the condition A £ H A is first satisfied in Model 2 as the frequency is 
decreased. It must be also noted that the derivative d\f\ 2 /dx is greatest in both models 
at the base which, according to our earlier discussion, indicates that locally the strongest 
wave reflection occurs there. The results presented in Figure 2e demonstrate similar wave 
behavior as that shown in Figure le. There is, however, a new interesting feature that can 
be seen for waves of period equal to 1 day. The dip observed at x « 1.2 is related to the 
fact that there is a change in the local density (as well as in Alfven velocity) scale height 
and that this change affects locally wave reflection; therefore, the derivative d\f\ 2 /dx goes 
through another zero (no reflection at this point). Figure 2f shows more convincingly than 
Figure If that after constant increase of the 7Z value with the period, there are frequencies 
for which reflection has a minimum and then starts to increase again. 

The described type of behavior is called “a tunneling effect” [e.g., Velli et al. 1993; 
Lou and Rosner 1993] to emphasize the fact that a travelling wave of a certain frequency 
can appear beyond the region of strong reflection as a propagating wave with some portion 
of its energy being lost due to reflection. To understand the origin of this effect, we now 
compare the results presented in Figures 2e and 2f. It is seen that only the plot of the 
reflection coefficient clearly indicates the existence of the tunneling effect. The plot of 
the function |/| 2 shows instead that non-WKB effects (stronger reflection) become more 
prominent for longer wave periods. Note that a similar discrepancy is briefly mentioned in 
the previous paragraph. The discrepancy between the results obtained by using |/| and 
'll is primary caused by the fact that the latter is defined in terms of \f\* and \g\ at the 
model base, and that values of these functions at the base are determined by the global 
wave behavior in the lower part of the model. The fact that this global behavior may be 
complicated is clearly shown in Figure 2e (see also Fig. le). Again, it is helpful to consider 
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the derivative of |/| 2 with respect to x, take its absolute value, and find the location of 
maxima and zeros. It is found that d\f\ 2 /dx for waves with periods of the order of 0.1 
day and shorter has only one maximum at the base. However, for waves with periods of 
the order of 1 day, there are clearly two maxima corresponding to regions of locally strong 
wave reflection; note that these regions are located symmetrically on both sides of the 
point where no reflection occurs ( dV A /dx = 0). If waves of certain periods are trapped 
between these two regions, then, changes in the global distribution of |/| (i.e, formation 

of a local maximum) are expected. As a result of those changes, the values of |/| and | <j| 
at the base become lower for long period waves than for short period waves (see Fig. 2e). 
This causes the wave reflection as measured by the value of the reflection coefficient to be 
stronger for short period waves. It is evident from the above discussion that this variation 
of the reflection coefficient with period is potentially misleading because in reality long 
period waves are more non-WKB waves than short period waves. This also means that 
the result may lead to a misinterpretation of Alfven wave reflection, in particular, if the 
reflection coefficient is solely used to study the wave behavior. 

After studying the influence of the magnetic field and density on the reflection of 
Alfven waves, we now proceed to the overall model dependence. We constructed four 
distinct wind models characterized by the same magnetic field at the base ( B 0 = 10G) but 
with the base temperature ranging from 7.45 x 10 5 K to 1.51 x 10* K (Models 4 - 7, see Fig. 
3a); this corresponds to the base density ranging from 5.0 x 10 15 <rm~* to 6.55 x 10 7 cm 
(see Fig. 3b). Although the high density at the base in Model 4 seems to be unrealistic for 
the Sun, we include this model into our parametric studies to see the effect caused by this 
extreme case on Alfven wave reflection (see Fig. 3b, for comparison to Sittler s model). 
As shown in Figures 3d and 3f, the efficiency of Alfven wave reflection in this particular 
model is very low but increases in wind models with higher temperature and lower density. 
Figure 3e presents similar behavior of the wave amplitude |/| as that observed in Figure 
2e. Waves with periods 0.04 d have their local maximum where waves with periods 5 d have 
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their local minimum. This is strong evidence for a complicated behavior of Alfven waves 
in both discussed models. In the light of our discussion of the results presented in Figures 
1 and 2, it is now easy to understand the behavior of Alfven waves of different penods in 
different solar wind models. The main conclusion is that Alfven waves are more efficiently 
reflected in models with the base temperature higher than 10* K and with the base density 
lower than 7 x 10 7 cm -3 . The results presented in Figure 3f also indicate that Alfven waves 
must have periods of at least several hours in cooler and denser atmospheres in order to 
be strongly reflected (71 > 50%). This is in contrast with the results obtained by H&O 
who showed that only waves with periods of one day or longer are strongly reflected in the 
solar wind. The difference can be easily explained by the fact that the efficiency of Alfven 
wave reflection is strongly model dependent. As we said earlier, using the global reflection 
coefficient to infer how freely waves propagate can be misleading. This is due to the fact 
(as explained in the previous paragraph) that the value of |/| at the base may become 
smaller even when periods increase. This again implies that the results based solely on the 
global reflection coefficient must be taken with caution. 

In all calculations described here, the condition \g\ 2 - |/| = const has been checked 
for each radial distance. We have found that the condition has always been satisfied with 
accuracy better than 1% in the inner (below the Alfvenic critical point) part of the solar 
wind model and that the accuracy is slightly lower far beyond the Alfvenic critical point. 

The discrepancy is due to the nature of Eqs. (5a, b), which are stiff for large values of a. 
However, this numerical check on our accuracy gives us full confidence in our results. 

WAVE ACCELERATION 

One of the main conclusions from the H&O paper is that non-WKB Alfven waves are 
no more effective in accelerating the solar wind than WKB waves. As discussed earlier 
in this paper, some recent results seem to lead to the opposite conclusion (Barkhudarov 
1991, Velli 1993). Here, we verify the H&O results by calculating the wave acceleration 
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for WKB and non-WKB Alfven waves. Following the formalism of H&O in deriving the 
formula for the wave acceleration, the proper way to calculate the acceleration force, which 
is caused by the propagating waves and acts on the background medium, is to consider 
the conservation of wave energy [e.g., Dewar 1970, Jacques 1977, H&O]. The wave energy 
conservation equation is 

^ + V«E= -fw-V , (12) 

at 

with the wave energy density, e, being defined by 

< = + 0 . (13) 
the wave energy flux, E, being given by 

E = pV^; + ^-r/ 1/ V^) * (14) 

and with the work done by the waves on the background medium being described by 

f.v=^{v; 2 -0-v-v|T (is) 


After taking the time average of Eq. (12) and using Eq. (14), we obtain 

f. . v = _ 1 • V . V f i v/ + v M 2 - Va*) . 

2 * 

which gives the expression for the wave acceleration 



(16) 


(17) 


where the averaged wave energy flux, E , is defined by 


E = l - P v l - V* + |Va *| 2 - 


Note that our formula for the wave acceleration (Eq. 17) differs from that used by H&O 
by a factor of 1/2 (see their Eq. 95). In our studies, Eq. (17) has been used to calculate 
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the wave acceleration in the wind models described in the previous section. The results 
are presented in Figure 4. We shall now discuss these results in detail. 

Our intention here is to investigate the dependence of wave acceleration on the mag- 
netic field strength at the base (Fig. 4a) and the wave period (Fig. 4b). The results 
presented in these figures clearly indicate that the wave acceleration decreases (i) when 
the magnetic field strength at the base increases and (ii) when Alfven waves with periods of 
several hours or longer are considered. As shown in the previous section, however, these are 
the cases when wave reflection increases and non-WKB effects become important. There- 
fore, we may conclude that the wave acceleration decreases for non-WKB waves. This 
rather surprising result can be explained by the following simple example. Let us consider 
a short period (say, several minutes) and a long period (say, several hours) Alfven wave 
entering the wind model at the base and carrying the same amount of energy; in the ap- 
proach used here the waves of different frequencies carry approximately the same energy. 
In all our wind models, the short period wave will propagate without any reflection but will 
exert the acceleration force on the background medium; the force is a local quantity and 
is proportional to the gradient of wave energy flux divided by the mass flux (see Eq. 17). 
The long period wave will undergo strong reflection and at a given height will always have 
much less energy available for exerting the acceleration force than the short period wave. 
This means that among WKB and non-WKB waves which carry comparable amount of 
energy, the former are always more effective in exerting a force on the background medium. 
The situation may be completely different, however, if the long period wave has a larger 
amplitude than the short period wave and thereby carries much more energy; the complete 
discussion of this problem is beyond the scope of this paper because it would involve a 
nonlinear wave treatment. 

The results presented and discussed above show why WKB waves are more effective 
in accelerating the background medium than non-WKB waves; it must be noted, however, 
that this result is valid only in a linear regime and when waves are initially assumed to 
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carry the same amount of energy. The problem we discuss now is the dependence of wave 
acceleration on different wind models. The results of our calculations are shown in Figures 
4c and 4d. In Fig. 4c, we plot the wave acceleration for the four different models described 
in the previous section. The main differences between the models are the values of physical 
parameters at the base (for example, in Model 4 the base temperature is the lowest and 
the base density is the highest; however, in Model 7 the base temperature is the highest 
and the base density is the lowest). In all models, we take the same magnetic field strength 
(B 0 = 2 G) and the same wave period (0.04* 1 ); it can be shown that Alfven waves with 
this period are WKB waves in all our models, since their function |/| 2 is either zero or 
stays close to zero everywhere in the model. Now, the results shown in Fig. 4c indicate 
that the wave acceleration is strongly model dependent (one could expect similar wave 
acceleration for the same WKB waves) and that it is the strongest for Model 7 and the 
lowest for Model 4. Because the results presented in Fig. 3f demonstrate that the efficiency 
of wave reflection increases from Model 7 to Model 4, with the exception of Model 6, one 
could therefore conclude that there is a contradiction between the results presented in Fig. 
4c (stronger reflection leads to stronger acceleration) and those shown in Fig. 4a and 4b 
(stronger reflection gives lower acceleration). To show that this conclusion is not correct, 
we perform the following calculations. We set the magnetic field strength and the wave 
period in such that the considered Alfven waves are primarily non* WKB waves in Model 
7 and almost WKB waves in Model 4; it must be remembered that by taking stronger 
magnetic fields we are allowing for more non- WKB effects. The results of our calculations 
are presented in Fig. 4d. First, we compare the wave acceleration in Model 4 and 7 
obtained for B a = 2 (?, and find a similar trend as that shown in Fig. 2c for Models 2 
and 3. Then, we compare the wave acceleration in Model 4 and 7 but for B 0 — 10 G to 
see that it has been reduced significantly by non- WKB effects and only slightly exceeds 
that obtained for waves in Model 4 that are nearly WKB waves. Finally, we compare the 
wave acceleration calculated in Model 4 for B a = 2 G (almost WKB waves) with that 
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obtained in Model 7 for B a = 10 G (strongly non-WKB waves). The comparison in Fig. 
4d shows that non-WKB effects reduce the wave acceleration and that it becomes even 
less effective than in Model 4; note that according to Fig. 4c, Model 4 has the lowest 
wave acceleration. In summary, the wave acceleration is strongly model dependent and it 
decreases with increasing non-WKB effects. 

The problem we want to address now is whether WKB waves are more effective in 
accelerating the solar wind than non-WKB waves. Despite the results we have obtained 
the problem is not trivial, because there is one additional condition that must be obeyed. 
The condition requires that the maximum of the wave acceleration occurs in the solar 
wind above the sonic critical point. As shown by Leer and Holzer (1980), only in this case 
can the deposited wave momentum be used to drive the solar wind to the observed large 
asymptotic flow speed. We have also studied this problem and the results are presented 
in Fig. 5, where we plot the location of the maximum of wave acceleration with respect to 
the location of the sonic critical point for two different wind models. The figure shows that 
longer period Alfven waves (primarily non-WKB waves) have maximum wave acceleration 
located farther above the sonic critical point than shorter period waves (primarily WKB 
waves). It is also seen that the range of wave periods which leads to this is extremely 
sensitively to the wind model. Our order of magnitude estimates (we use the H&O method 
of converting dimensionless variables into their physical values) indicate that in most cases 
even WKB waves exert a force which is too small to accelerate the wind; the obtained result 
is, however, affected by the constraint that the considered Alfven waves carry roughly the 
same amount of energy despite their periods. In addition, as we clearly presented in 
Figs. 3 and 4, all the results are strongly model dependent and it can be shown that for 
some cases the wave acceleration can exceed the wind acceleration before the sonic critical 
point. Therefore, one can expect that Alfven waves could play an important role in the 
atmospheres of other stars. 
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CONCLUDING REMARKS 


In this paper, we discussed the behavior of non-coupled, undamped, linear and toroidal 
Alfven waves. Our studies indicate that both wave reflection and wave acceleration are 
strongly model dependent so that it is difficult to predict the atmospheric response to the 
waves without an appropriate wind model. The range of models considered in this paper 
allows us to draw several important conclusions. They can be summarized as follows: 

(1) Alfven waves undergo substantial reflection, providing that base magnetic field strength 
B 0 is at least 10G and/or the wave period is in the range of hours or longer. 

(2) The strongest reflection can be expected in solar wind models with the base tempera- 
tures higher than 10® K and with the base densities lower than 7 x 10 T cm -3 (neither of 
these values contradict presently known observational data). In these models as well as 
in the models with lower temperatures and higher densities, Alfven waves with periods as 
short as several minutes can be treated as WKB waves; for Alfven waves with periods of 
the order of one hour or longer a non- WKB treatment is required. 

(3) The region of strongest wave reflection is usually located at the base of the models. 

(4) The variation of the reflection coefficient with period can be misleading because it often 
shows that short period waves are more reflected than long period waves; in reality the 
latter are always more non- WKB waves than the former. 

(5) Finally, non- WKB, linear Alfven waves are always less effective in accelerating the 
plasma than WKB Alfven waves. 

In addition to these conclusions, we have attempted to give the physical reasons for our re- 
sult. In the following, we shall compare the presented results to those previously obtained. 

We begin with the comparison to the H&O results. The main difference is that 
our studies demonstrate that both Alfven wave reflection and Alfven wave acceleration 
are strongly model dependent. In addition, we show that in some wind models Alfven 
waves with periods as short as one hour must be treated as non- WKB waves. Our results 
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concerning the wave acceleration are in agreement with those obtained by H&O but they 
are in contradiction with some results obtained by Barkhudarov [1991] and Velli [1993] 
who found that non-WKB Alfven waves of certain frequencies can be more effective in 
accelerating the plasma than WKB waves; this is exactly opposite to one of the main 
conclusions of this paper. Velli [1993] also calculated the dependence of the transmission 
coefficient of Alfven waves on their periods and found that it is very sensitive to the 
temperature both in a static atmosphere and in an isothermal solar wind. In the case of 
a static atmosphere, only waves with periods less than about 15 minutes are completely 
transmitted whereas in the case of the moving atmosphere, he showed that the transmission 
coefficient generally increased even for periods in the range of tens of days. He claimed that 
this type of behavior is due to the changes in the temperature in his isothermal models. 
We also differ from Velli in that we show significant reflection for long period waves; his 
transmission coefficient almost equals unity for these periods. This is likely explained by 
using his Eq. (25) and the adjustment of the value of the Alfven velocity at the solar 
surface and at the Alfvenic critical point to nearly the same value, which is not the case 
in our studies. This tuning of Alfven velocities might be an artificial effect of computing 
his Alfven velocity profile and may also affect his calculations of the wave acceleration. In 
general, Velli and his co-workers have drawn a number of conclusions concerning Alfven 
wave reflection by using the transmission (reflection) coefficient. As shown in this paper, 
interpretation of results based solely on the concept of the reflection coefficient can be 
misleading because it often reveals that short period waves are more reflected than long 
period waves; in reality the opposite is true. 

Finally, our results related to the effectiveness of reflection for short period waves 
seems to be in contradiction with the previous claims made by Moore et al. [1991]. They 
found that for a static, isothermal and radially spreading coronal hole with 2?® — 10 G 
at r 0 = 1.15 Rq all waves with periods longer than several minutes should be reflected if 
the temperature at the base is less than 1.0 x 10® /if, and should escape to the solar wind 
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if that temperature is higher than 1.0 x 10® K. They suggested that these short period 
(reflected) Alfven waves are responsible for the heating in solar coronal holes. The results 
obtained in this paper for non-isothermal models clearly indicate that Alfven waves with 
periods as short as several minutes cannot be reflected in the solar wind models; similar 
results have been obtained when we considered isothermal models. Therefore, we conclude 
that the approach presented by these authors must be revised. 
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FIGURE CAPTIONS 


Fig. 1. Reflection of Alfven waves in the solar wind Model 1 (polytropic index a p 1.05 
here and in all succeeding models), a) Temperature distribution with T a = 1.63 x 10 «K. b) 
Electron density (n eo = 8 x 10 7 cm" 1 ) compared with that used by Heinemann & Olbert, 
called “Sittler’s model”, c) Different values of the magnetic field at the base of corona, 
2, 5, 10 and 15 Gauss, respectively, result in different profiles of Alfven velocity. Solid 
line denotes a wind velocity consistent with the calculated atmosphere model, d) Radial 
dependence of the wave amplitude |/| 2 for magnetic fields corresponding to these in panel 
(c). e) A case study of |/| 2 versus distance for waves of different period with a constant 
value of B 0 = 10G at the base of the model, f) The reflection coefficient, TZ, as a function 
of wave period for four different magnetic field strengths, B 0 . 


Fig. 2. Reflection of Alfven waves in the solar wind Models 2 and 3, with differing base 
density but nearly the same temperature, a) The temperature profiles T 0 = 1.28 x 10 
and T 0 = 1.35 x 10® for Model 2 and Model 3. respectively, b) The density profiles. Initial 
electron densities are 7.0 x 10® and 3.75 x 10® cm" 3 , respectively, c) The wind and Alfven 
velocities for both models, d) The wave amplitude i /| 2 in these models for field B a = 10 
G and period r = 0.04*. e) Wave reflection versus distance for Model 2 for different wave 
periods, f) The reflection coefficient, 71, is plotted as a function of wave period for both 

considered models. 


Fig. 3. Alfven wind reflection in wind models with strongly differing base density and 
temperature, a) Temperature distributions for the four models with the extreme ones 
differing by a factor of two (Model 4 and Model 7, respectively), b) Density profiles. The 
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difference between respective electron densities spans 5 orders of magnitude, with higher 
temperature corresponding to the lower density, c) Alfven speeds and flow speeds. Despite 
the huge differences between densities, Alfvenic critical point locations and corresponding 
critical velocities, are concentrated around similar values. The magnetic field is 10 Gauss 
for all cases, d) The wave amplitude |/| 2 for the four models, e) Dependence of |/| 2 on 
the wave period in Model 5 with a fixed magnetic field of 10 Gauss at the base of the 
atmosphere. The maximum for wave with r = 0.04** and local minimum for wave with 
r _ 5<i are located at the point where density distribution changes its dependence with 
the distance, f) The reflection coefficient, K, as a function of period for the four models. 


Fig. 4. The dimensionless wave acceleration as a function of radial distance from the Sun 
for several models, a) Effect of varying field strength in Model 7 (see Fig. 3). b) Effect 
of varying wave period in Model 7. c) Models 4 - 7 compared, d) Models 4 and 7 for the 
field strengths 2 Gauss and 10 Gauss. 


Fig. 5. The difference between the location of the point where acceleration has its maxi- 
mum and the location of the sonic point (where the wind velocity equals to the local sound 
speed) as a function of wave period. This difference is plotted for Models 5 and 7 (see Fig. 

3). 
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